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We investigate the dephasing rate, 1/r^,, of weakly disordered electrons due to scattering from diluted 
dynamical impurities. Our previous result for the weak-localization dephasing rate is generalized from 
diluted Kondo impurities to arbitrary dynamical defects with typical energy transfer larger than l/r v . For 
magnetic impurities, we study the influence of magnetic fields on the dephasing of Aharonov-Bohm oscil- 
lations and universal conductance fluctuations both analytically and using the numerical renormalization 
group. These results are compared to recent experiments. 
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I. INTRODUCTION 

Decoherence is the fundamental process leading to a sup- 
pression of quantum mechanical interference and therefore 
is indispensable for our understanding of the appearance of 
the classical world. The destruction of phase coherence in 
a quantum system occurs due to interactions with its en- 
vironment and can be studied, e.g., in mesoscopic metals 
and semiconductors where the quantum-mechanical wave 
nature of the electrons leads to a variety of novel transport 
phenomena at low temperatures. 

Although the concrete definition of the dephasing rate, 
1/t^, depends on the experiment used to determine it, the 
electron-electron interactions are thought to be the dom- 
inant mechanism for the destruction of phase coherence 
in metals without dynamical impurities below about 1 K. 
The dephasing rate for interacting electrons in a diffusive 
environment was first calculated by Altshuler, Aronov and 
Khmelnitsky (AAK) and vanishes at low temperatures, T, 
with some power of T, depending on the dimensionality of 
the system^ 

In the last decade several independent groups performed 
magnetoresistance experiment s 2 ' 3 ' 4 ' 5 to probe the influence 
of dephasing on weak localization in disordered metallic 
wires. Irritatingly, a saturation of the dephasing rate, 
has been observed at the lowest experimentally accessible 
temperatures. This observation has triggered an intense 
discussion on the mechanism responsible for the excess of 
dephasing. 6,7,8 The most promising candidates to explain 
the saturation of 1/t v are extremely low concentrations of 
dynamical impurities, such as atomic two-level systems 9,10 
or magnetic impuritie s 3 ' 4 ' 5 ' 11 ' 12 ' 13 ' 14 . This has been corrobo- 
rated on the one hand by experiments 3,5 on extremely clean 
Ag and Au samples where the dephasing rate continues to 
decrease well below 100 mK and on the other hand by dop- 
ing studies with magnetic impurities. 2 i 3 i 5 i 15 i 16 As expected 
theoretically, 12 these experiments show a saturation of 1/t v 
above the Kondo temperature, Tk, the characteristic scale 
of screening of the magnetic moment, and a suppression 
of I/Tip below this scale. Recent highly controlled exper- 
iments,^ 5 -^ in which a few ppm (parts per million) of Fe 
ions have been implanted by ion beam lithography into very 
clean Ag samples, showed that the screening of these Fe 



ions is surprisingly well described 15 by the theoretically pre- 
dicted dephasing rate for spin-1/2 Kondo impurities 13 down 
to temperatures of 0.1 Tk. At the lowest temperature, again 
a plateau in the dephasing rate has been observed propor- 
tional to the number of implanted Fe ions. The origin of this 
puzzling behavior is still unclear but may arise from further 
dynamical defects created during the implantation process 
or by rare Fe ions with a different chemical environment and 
strongly reduced magnetic screening. 

An obvious option to study the influence of magnetic 
impurities on the dephasing rate is to measure its depen- 
dence on an externally applied magnetic field. The applica- 
tion of sufficiently large magnetic fields freezes out inelastic 
spin-flip processes as discussed theoretically in Refs. fl2j ]. 
and therefore one expects the dephasing rate to return 
to the value predicted by AAK for dephasing induced by 
Coulomb interactions in a diffusive environment. The orbital 
contribution of the magnetic field does, however, destroy 
the weak-localization (WL) contribution to the magneto- 
resistance, as the joint propagation of an electron and a 
hole along time-reversed trajectories (the Cooperon) picks 
up extra (random) Aharonov-Bohm phases in the presence 
of external magnetic fields. Measuring the independent de- 
phasing rate in a WL experiment is therefore only possible 
in strictly one- or two-dimensional systems using magnetic 
fields almost exactly parallel to such a structure, requiring 
an accurate alignment of magnetic fields. 

Universal conductance fluctuations (UCF) and Aharonov- 
Bohm (AB) oscillations with a periodicity of h/e, on the 
other hand, are not suppressed by orbital effects and can 
be used rather directly to determine the field dependence of 
the dephasing rate. 17 

UCFs can be observed as characteristic fluctuations of the 
conductance as a function of the magnetic field. The exter- 
nal magnetic field enters the metal and changes the pattern 
of the electrons wave functions and therefore the conduc- 
tance in a random but reproducible way ( "magnetofinger- 
print"). In AB experiments performed on mesoscopic rings, 
these sample fluctuations are further modulated by periodic 
/i/e-oscillations resulting from the magnetic flux piercing 
the ring. Both UCF and AB oscillations rely on the con- 
structive interference occurring in the collective propaga- 
tion of electrons and holes traveling along the same path 



(the diffuson). These are robust against the breaking of 
time-reversal invariance (while Cooperon contributions are 
rapidly suppressed by small fields) but are sensitive to de- 
phasing by inelastic processes. Indeed, Benoit et a/3 and 
Pierre et a/. 5 have shown that the amplitude of Aharonov- 
Bohm oscillations increases by almost an order of magnitude 
for increasing magnetic fields clearly showing a suppression 
of dephasing by magnetic (Zeeman) fields. This leads to 
the conclusion that the main mechanism of dephasing in 
the investigated low-temperature regime is the scattering 
from magnetic impurities. 

Previously we have studied the zero-field dephasing rate 
due to diluted Kondo impurities as measured from the WL 
experiment. 13 We showed that the dephasing rate for all ex- 
perimentally relevant temperatures is proportional to the in- 
elastic crosssection, 14 which itself can be expressed in terms 
of the T-matrix describing the scattering of the electrons 
from a single magnetic impurity. Such a relation has been 
proposed previously by Schwab and Eckern 19 in the context 
of UCFs. 

In this paper we generalize our previous work to arbitrary, 
diluted impurity scatterers with typical energy transfer larger 
than l/r„. Furthermore we supply the magnetic field de- 
pendence of 1/t„ as measured in the AB experiment (and 
compare our results to the AB experiments performed by 
Pierre and Birge. 20 ) The outline of this paper is as follows: 
In Sec. II we discuss the dephasing rate due to generic dy- 
namical scatterers as measured in the WL experiment. We 
briefly review our previous results 13 and generalize them to 
arbitrary dynamical impurities with typical energy transfer 
larger than In Sec. Ill we turn to the dephasing rate 

as measured from the UCF and the amplitude of the AB os- 
cillations. We briefly review the main concepts entering the 
analysis of these experiments and discuss how the dephasing 
rate measured from the UCF differs from that measured in 
the WL experiment. The main goal of Sec. Ill is to give the 
dephasing rate as measured from the amplitude of the AB 
oscillations. Results for the dephasing rate obtained using 
the numerical renormalization group (NRG) are described 
in Sec. IV. Sec. V summarizes with a discussion. 



see Fig. [TJ In the absence of dephasing by inelastic pro- 
cesses Co(q) is the bare Cooperon, 
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Dcf + in + 1/tb ' 
as diagrammatically depicted in Fig. [2]and 21 
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is the dephasing rate due to the applied magnetic field, B. 
Here D is the diffusion constant, d is the dimension of the 
diffusion process, r denotes the mean scattering time cor- 
responding to a mean free path / = vft, n is a unit vector 
orthogonal to the probe in d = 2 and pointing along the wire 
in d = 1 and L± is the transverse dimension of the sample. 
As can be seen from Eq. JT} the WL corrections depend 
on the strength of the applied magnetic field, B, and di- 
verge in low dimensions, d = 1, 2 for B = 0, reflecting the 
fact that WL corrections in low-dimensional systems may 
become strong and lead to strong Anderson localization. 

Taking into account interactions (as e.g. provided by dy- 
namical impurities) the bare Cooperon dresses with a mass, 
1/t„, i.e. (if purely exponential decay is guaranteed), 



Cn =0 (q) 



Dqi + 1/t v + 1/t b ' 



(5) 



For weak magnetic fields (t v <C tb) the WL corrections 
are therefore cut off by 1/t v , allowing to determine the 
dephasing rate from fitting the magnetoresistance. 

In this section we study the dephasing rate due to generic 
diluted, dynamical scatterers with typical energy transfer 
larger than 1/t v as measured from WL. To be specific, we 
consider a Hamiltonian of the general form 



^imp — -^imp + c kcr C k'o-'/kk'(T(T'^o 



i(k-k')x; 



(6) 



II. DEPHASING RATE FROM WEAK-LOCALIZATION 
CORRECTIONS TO THE CONDUCTIVITY 



The most accurate way to extract l/r„ at low mag- 
netic fields is via the WL corrections to the Drude con- 
ductivity, which result from coherent back scattering of an 
electron-hole pair traveling along time-reversed paths in the 
disordered environment. Technically, the coherent propaga- 
tion of the electron-hole pair is described by the Cooperon, 
Co(q), and the WL correction is given by 



2e 2 D 



(1) 



where Hf m? is the Hamiltonian of the isolated impurity, c', c 



are creation and annihilation operators of conduction band 
electrons, Xj denotes the position of the impurities, and the 




FIG. 1: Diagrammatic representation of the Cooperon 
Cn =0 (q) which enters the WL corrections to the Drude con- 
ductivity. p± = q/2 ±p and wavy lines denote current oper- 
ators. 
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FIG. 2: Bethe-Salpeter equation for the bare Cooperon 
Co(q, SI) = 2^C (q, SI), p± = g/2 ± p, and p'± = q/2 ± p'. 
Dashed lines denote scattering from static impurities and R, 
A denotes the particle and hole lines, respectively. 



momentum and spin dependent function f a parametrizes 
the coupling to some operator X a describing transitions of 
the internal states of the dynamical impurity. Eq. ([6]), e.g., 
describes the coupling of the conduction band electrons to 
Kondo impurities, two-level systems, etc. In order to find 
the dephasing rate due to such generic impurities one has 
to compute the "self energy" or "mass" of the Cooperon 
generated by the operators of Eq. ((6|). 

Two assumptions allow to reduce this problem to that of 
summing up a simple geometric series. First, we assume 
that the concentration of dynamical impurities rii is small, 
and, second, that the elastic mean free path I is large. How 
small rii and how large I has to be, depends on both, the 
dynamics and the extension of the impurity as briefly dis- 
cussed below and - for Kondo impurities - in more detail 
in Ref.[13] (the influence of stronger disorder on the dis- 
tribution of Kondo temperatures has recently been investi- 
gated by Kettemann and Mucciolc 22 ). Note that in relevant 
experimental systems in the WL regime 2,3,4,5 the two as- 
sumptions are well justified^ 3 -. The first observation is that 
quantum interference corrections to the inelastic scattering 
rate are small when a diffusing electron is unlikely to return 
to the same dynamical impurity (i.e., when the weak local- 
ization corrections are weak). Technically speaking, this is 
reflected in the fact that diagrams mixing scattering from 
dynamical and static impurities are suppressed by factors 
of 1/N± and a/l [or l/(kpl) for a < where N± is 

the number of transverse channels in a quasi 1 or 2 dimen- 
sional system and a is the typical diameter of the dynamical 
impurity. This effect and further system dependent factors 
relevant for the suppression of quantum interference correc- 
tions are discussed in Ref.(l3|. Only at lowest, experimen- 
tally unprobed temperatures, does the enhanced infrared 
singularity, caused by the presence of extra diffusion modes, 
overcompensate this phase space suppression factor, as dis- 
cussed in detail in Ref. [l3j]. The small parameter l/(kpl) or 
a/l therefore reduces the problem to compute the 'mass' of 
the Cooperon to that of solving the Bethe-Salpeter equation 
diagrammatically depicted in Fig. [2a). 

For small rii, one can furthermore restrict the analysis of 
the irreducible vertex T to terms linear in rii as shown in 
Fig- [3jb). T can be separated into three distinct contribu- 
tions: self-energy diagrams [the first two terms in Fig.[3{b)], 
an 'elastic' vertex correction with no energy transfer be- 
tween upper and lower line (third term), and an 'inelastic' 
vertex where interaction lines connect the two lines (last 
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FIG. 3: Bethe-Salpeter equation for the Cooperon C in the 
presence of dilute dynamical impurities to linear order in m. 
Co is the bare Cooperon in the absence of interactions and T 
the irreducible vertex obtained by adding self-energy, elastic-, 
and inelastic- vertex contributions. The crosses with attached 
dashed lines denote the averaging over impurity positions Xi, 
the squares the inelastic scattering from a single impurity to 
arbitrary order. 



term). Only this inelastic vertex makes the Bethe-Salpeter 
equation a true integral equation as it mixes frequencies 
but, fortunately, this term can be neglecteo- if the typical 
energy, AE, exchanged between electrons and holes during 
an interaction process greatly exceeds the dephasing rate 
due to the dynamical impurities, l/r v , i.e., AEt v 3> 1. 
Physically, 1 the suppression of the inelastic vertex arises as 
an exchange of energy AE leads to a phase mis-match 
of order e lAETv between electron and hole destroying in- 
terference completely for AEt v ^> 1. Technically, one 
can confirm this argument by estimating corrections to the 
WL contributions due to the inelastic vertex, as, e.g., de- 
picted in Fig. [4j 13 More importantly, however, the condition 
AEt v 3> 1 always holds for sufficiently small concentra- 
tions rii, since 1/t v cx rij. In the case of Kondo impurities 
discussed below, for example, this condition translates to 
m < vT K . 

Restricting to the self-energy and elastic vertex contribu- 
tions, the Bethe-Salpeter equation is easily solved: Since 
self-energy and elastic vertex contributions conserve the en- 
ergy of single electron lines, the solution of the reduced 
Bethe-Salpeter equation amounts to a straightforward sum- 
mation of a geometric series. Setting the center-of-mass 
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FIG. 4: Lowest order correction to WL contributions due to 
inelastic vertex. The Cooperon, C, is dressed with a mass 
resulting from summation of the elastic part of the vertex F, 
i.e., the self energy and the elastic vertex contribution. 



frequency ft to (see Fig. [JJ and Fig. [3]), the Cooperon is 
given by 



C n =o(e,q) = 



£q 2 + l/^(e, T) + 1/t b ' 
with the T and e dependent dephasing rate 
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restricts the electrons mo- 



Here ,g e (p) 

menta, p, and energies, e, to the Fermi-surface, e(p) is 
the dispersion relation of the conduction band, T A,R are the 
advanced /retarded T-matrices, defined by the Green func- 
tion G xx ,(e) - G° x ,(e) + G x0 (e)T(e)G° x , (e), and v de- 
notes the density of states per spin. Eq. (jS} generalizes our 
result of Ref. [1J| to arbitrary shaped diluted impurities. 
Notice that also forward scattering processes enter 1/t„, 
which do not contribute to the transport scattering rate. 
We stress that Eq. ([8]) is the general result for the dephasing 
rate for a weakly disordered metal due to a low concentra- 
tion of generic dynamical impurities for which the condition 
AE > \/t<p holds. In the opposite limit, AE > 1/t v , 
vertex corrections become important, as has been discussed 
earlier on 23 ' 24 in the context of magnetic impurities. 

As we assumed that a < I, Eq. ([8]) can be further sim- 
plified, 
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where Sp is the Fermi-surface (or more precisely the surface 
with efc = e). Here we also assumed a time-reversal invari- 
ant system with TJL(e) = T R p , p (e) and employed the 

identity [Tp* ,(e)]* = T$ (e). Also Kettemann and Muc- 
ciolo have generalized the dephasing rate for a momentum 
independent T matrix to Eq. ([9]) independently in a recent 
report. 22 The dephasing rate given in Eq. ([9]) has a simple 
interpretation-! 4 -: Since the Fermi-surface integrated imagi- 
nary part of the T-matrix is proportional to the total cross- 
section and |T pp /| 2 (integrated over the Fermi surface) is 
proportional to the elastic cross section, its difference is, by 
definition, proportional to the inelastic cross section, <7i nc i, 
introduced in Ref. Therefore, Eq. © can be rewritten 
in the form 



1 



7V, (e) 



= ni(v F (p) CT inc i(p, e)) , 



(10) 



where (...) denotes an angular average weighted by l/vp(p) 
to take into account that fast electrons are scattered more 
frequently from elastic impurities. According to Eq. (|10p . 
t v is nothing but the average time needed (in a semiclassi- 
cal picture) to scatter from an impurity with cross section 
CTinci- Note that the vanishing of 1/t v for static impurities 
is guaranteed by the optical theorem. 

From Eq. ([9]) we can read off the dephasing rate for di- 
luted dynamical isotropic s-wave scatterers 13 19 
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where v loc is the local density of states at the Fermi en- 
ergy at the site of the impurity which can differ from 
the thermodynamic density of states entering the prefac- 
tor. Note that in the case of Kondo impurities discussed 
below (and in Ref. 0), the combination v Xoc T R / A {t) = 
/(e/Tk, T/Tk, B/Tk) is an un/Versa/dimensionless function 
of the ratios e/Ty^,T/Ty^, B/T^. If the assumptions under- 
lying the derivation of Eq. (jlip are valid, one can therefore 
predict without any free parameter the dephasing rate if the 
concentration of spin-1/2 impurities, the Kondo tempera- 
ture and the thermodynamic density of states are known 
(see, e.g., Ref. [15]). However, one of the assumptions un- 
derlying the derivation of the prefactor of Eq. (fTTjl may not 
be valid in realistic materials: we assumed that the static im- 
purities are completely uncorrelated and local such that elec- 
trons are scattered uniformly over the Fermi surface. While 
this should be a good assumption in doped semiconductors, 
this may not be valid in metals with complex Fermi surfaces 
and strongly varying Fermi velocities. Under the latter con- 
ditions, we expect that the prefactor of Eq. (fTTj) becomes 
nonuniversal, yielding temperature-independent corrections 
of order one, which may be important for the interpretation 
of high-precision experiments. 15 ' 16 

In Ref. [IH we have calculated the leading corrections 
to Eq. (fTTj) arising from mixed diagrams involving com- 
bined scattering from static and dynamical impurities and 
from diagrams including higher processes in rij. We showed 
that, suppressed by the small parameter l/(kfl), their con- 
tributions are negligible at all experimentally relevant tem- 
peratures, T. Only at the lowest experimentally irrelevant 
temperatures do these corrections become important due 
to infrared singularities of the dressed interaction potential 
(dressed by coherent backscattering processes). The esti- 
mates of subleading corrections presented in Ref. [l3[ can 
be generalized to extended dynamical impurities by replac- 
ing l/(kpl) by a/l for kpa > 1. 

The experimentally measured dephasing rate does not re- 
solve the dependence on the electrons energy e. We de- 
scribed in Ref. [13] that to allow for a comparison with 
the e-independent dephasing rate, r ¥ 7 1 (T), extracted from 
the WL experiment, the energy-resolved representations of 
1/tu,, Eq. ([9])- (fTTjl. still require an average over energies 
according to 



Jdef^)r v (e,T)^ 

7de^(e)ln^ 
fdef>(e)/r v (e,T) 



^exp 



d= 1,3, 

d = 2, 
t>/t s > 1. 



(12) 



Here the last line applies to a case where a relatively strong 
magnetic field, B, is present. 

Specifying to a situation where the coupling of the con- 
duction band electrons to the (diluted) dynamical impuri- 
ties is described by the spin-1/2 Kondo effect, the general 
Hamiltonian of Eq. ([6]) takes the form 



Hs = J Sjcl (Xj )a rTCT ' c a > (Xj) , 



(13) 



where J is the exchange constant. An external magnetic 
field, B, causes Zeeman splitting, e z , of the conduction 
band electron spin states, 



e 2 = 9c^bB, 
and couples to the impurity spins according to 



(14) 



For wires of length L Lt = \J D/T, the fluctuations 
of the conductance, 5g5g, are determined 25 by 



5 ^9=^^ / </<^,/F(n>/F(-) 



J dxidx 2 |P ei ,e 2 (xi,X2) 
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Here P £lj£2 (xi,X2) is the amplitude for an electron-hole 
pair, with energies e%,62 respectively, to diffusively travel 
from xi to X2 along the same trajectory (diffuson, see 
Fig. EI). The overbar denotes the ensemble average, which 
is experimentally realized by changing the magnetic field. 
Eq. (|16p assumes the large ring diameters, L ^> L v , such 
that the dephasing rate, l/r^, controls the magnitude of 
the fluctuations. Furthermore we assume temperatures 
Tt v ^S> 1. Notice that generally there is also a contribution 
from the Cooperon, which is, however, suppressed already 
for small magnetic fields. 

We change to momentum representation and separate the 
two-particle propagator P into its spin-singlet and -triplet 
components, P 
tation of Ref. [l 



Yli—i 4-P i where (following the no- 
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g e , gs are the electrons and the magnetic impurities gyro- 
magnetic factors, respectively. As already mentioned in the 
introduction, measuring the £?-field dependence of the de- 
phasing rate due to Kondo impurities in a WL experiment 
is a highly delicate task. In view of this difficulty it is more 
feasible to measure the independence of 1/t v from the am- 
plitude of the AB oscillations as discussed in the following 
section. 



I. DEPHASING RATE FROM UNIVERSAL 
CONDUCTANCE FLUCTUATIONS AND 
AHARANOV-BOHM OSCILLATIONS 



(17) 



The modes i = 1,2,3 describe the spin triplet state with 
S z component equal to 1,-1, and 0, respectively, i = 4 
denotes the spin singlet channel. The Zeeman splitting en- 
ters only the triplet-diffuson with nonvanishing projection 
S z = ±1, i.e. (i = il f° r i = 1, 2 and zero otherwise. 
V T so ' s tne spin-orbit scattering rate. 1/tjq is identical 
for the three spin triplet-diffuson (i = 1,2,3) and zero for 
the spin singlet mode (i = 4). For strong spin-orbit scatter- 
ing only the singlet diffusion contributes (otherwise 1 / T so 
is an additional fitting parameter). Finally 1/t^ s is the de- 
phasing rate for the zth diffuson mode due to the presence 



Let us first begin with a brief discussion of the universal 
conductance fluctuations (UCF) and their dependence on 
1/t„ and then turn to the experiment on Aharonov-Bohm 
rings. 

To be specific, we consider a wire of non-interacting elec- 
trons, scattering elastically from static impurities, and in- 
elastically off a low concentration, rii, of Kondo impuri- 
ties, where the coupling of the conduction band electrons 
to the dynamical impurities is described by the Hamiltonian 
Eq. (fT3|) and the influence of the magnetic field is accounted 
for by Eqs. (HID and (JIHJ). 
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FIG. 5: Bethe-Salpeter 
D {q,n) = ^D (q,Cl). 
Dashed lines denote scattering from static impurities and R 
A denotes the particle and hole lines, respectively. 



equation for the bare diffuson 
p ± = p ± q/2 and p'± = p ± q/2. 



R 




FIG. 6: Diagram giving the main contribution to the UCF. 
Dashed lines represent coherent impurity scattering of elec- 
tron (R) hole (A) pair, i.e. the bare diffuson where interaction 
due to scattering from magnetic impurities is not yet taken 
into account. 



of diluted magnetic impurities which has the structure 

W-^ = ^fe [T(<,B)(ea * fi) - T(,,6)(ei '- B)] 

-(7r^) 2 T^( ei ,i3)T^(e 2 ,i3)^ (18) 

where the proper combination of T-matrices for the vari- 
ous channels can be read off by comparison with Table [T] 
Eq. (fl8|) is evaluated from summing up self-energy and elas- 
tic vertex contributions. Notice that in contrast to the 
WL experiment the electron and hole lines (i.e., the inner 
and outer rings) in Fig. [6] represent different measurements. 
Therefore there are no correlations between dynamical im- 
purities residing on different rings and interaction lines may 
only be drawn within the same ring. Consequently the in- 
elastic vertex contributions do not enter the Bethe-Salpeter 
equation for the diffuson, see Fig. [JJ Notice that there are 
inelastic vertex contributions, as, e.g., depicted in Fig. [8] 
which become important in the context of electron-electron 
interactions^ It is instructive to compare those to the in- 
elastic vertex corrections relevant for WL depicted in Fig. 
|4] In the latter case, the sum of the incoming momenta of 
the vertex is small due to the Cooperon in Fig. [4] Conse- 
quently, the inelastic vertex corrections to the WL dephas- 
ing rate are not suppressed by powers of l/(kfl) but only by 
powers of 1/(ATt v ). In contrast, the relevant momenta in 
Fig.[5]are uncorrelated (i.e., particle and hole are far apart), 
leading to an suppression both by powers of l/(kpl) and of 




FIG. 7: (a) Bethe-Salpeter equation for the diffuson, D, in 
the presence of (dilute) magnetic impurities to linear order 
in rii. Do is the bare diffuson in the absence of interactions, 
(b) Diagrammatic representation of the irreducible interac- 
tion vertex, F, consisting of the self-energy (represented by 
the first two contributions), the elastic vertex (third contri- 
bution). The inelastic vertex (the fourth contribution) does 
not enter the diffuson as measured in the UCF. 



We point out the following differences for 1/t v .s mea- 
sured from the UCF experiment, Eq. (|18p. compared to 
that found from the WL, Eq. ([TTj). First, the T-matrices 
entering Eq. (|18p depend on the spin configuration of the 
diffuson-mode and have acquired a independence due to 
the coupling of the impurity spin to B, Eq. (fl5|) . Second, 
1/ T v-S depends on two energies. This results from the fact, 
that in the UCF experiment electron and hole lines consti- 
tuting the diffuson are produced in different measurements 
of the conductance (see Fig. [6}. Therefore their energies 
are individually averaged as can be seen in Eq. (fl6|) . 

From Eq. (|16p and (fT7|) the amplitude of the UCFs is 
obtained to be proportional to y/r v _s- Especially compared 
to the AB oscillations, discussed below, the depen- 
dence of the UCFs is rather weak. In the following, we will 
therefore focus our discussion on AB experiments. 

Aharonov-Bohm oscillations are measured in a ring ge- 
ometry ) 20 ' 28 where the conductance oscillates periodically 
as a function of B piercing the ring, and the amplitude of 
these AB-oscillations decrease exponentially with l/r Vt s- 
The periodic oscillations result from the change of bound- 
ary conditions, due to flux lines piercing the ring and can 



I 


\S,M) 


Combinations of T-matrices 


1 


S=1,M = 1 


Ti = ±(T?-T«)-T?T* 


2 


S = 1,M = -1 


T 2 = i (T A - Jf) - TfTf 


3 


S=1,M = 


T 3 = ± Im (T T A + T A ) - i (T T R T T A - Tf T A ) 


4 


s = o 


T*= I Im (7? + Tf) - 1 (T«T*-T«T*) 




TABLE I: Combination of T-matrices entering the dephasing 
rates for spin-triplet and spin-singlet diffusons. S denotes 
the total spin and M its z component. T|,T denotes the 
T-matrix for spin-up and spin-down electrons, respectively. 



FIG. 8: Diagrammatic representation of lowest order correc- 
tions to UCF due to inelastic vertex contributions. Notice 
that for a local interaction these are not only suppressed by 
powers of l/(ATr v ) but also small in powers l/(fcpi). 



be calculated from 



J dx 1 dx 2 |P^ 2 (xi,x 2 )| : 



(19) 

Here is again given by Eq. (fT7|) but now the con- 

tinuous q have to be replaced 25 by discrete momenta, 
Q = Qm(A0) = ^-(m + depending on the differ- 

ence of the magnetic flux during the individual measure- 
ments of g. The fluctuations are a periodic function in 
A(j)/(j>Q, where (f>o = 27r/e is the elementary flux-quantum 
and A<fi = ABL 2 /(4tt). Therefore an expansion in its har- 
monics can be made, 25 



6g((f>)5g(<f> + A<j>) 



-^-> A S cos 

k=0 



2,k^ 



(20) 



where C is a factor of order 1, depending on the sample 
geometry in the vicinity of the ring. Restricting to the 
situation of strong spin-orbit scattering 29 where the spin 
singlet-diffuson gives the leading contributions to Eq. (|19p , 
one finds that (for L ^> L^) 



A(B) 



(2tt) 3 Z? 3 / 2 



T 2 L 3 



de 



cosh 4 (e/2T) 



where \t v = rj^ s in the notation of Eq. (|18p] 



v,s 



1 



r v (e,T,B) 



2rii 



■kv Im 



Tf 4) ( e ,B)l -|(7ri/)Tf 4) (e,fl)| 5 



(22) 



Here e = ei + e 2 and we used that relevant contributions to 
the integral over energy differences, e = t\ — e 2 , result from 
energies e <C 1/t^ to eliminate the e dependence. Notice 
that such a reduction to a single energy-integral can only be 
done in a one-dimensional system where the q-integral over 
the square of the diffuson, Eq. (fT7|) , is dominated by infrared 
divergences. In a 2-d system, e.g., relevant energies extend 
to e ~ T. For the following, it is convenient to rewrite 
Eq. {22} [see also Eq. {TO}] as 



1 1 (Cinel) 



7)1 it ^max 



(23) 



where cr max = An/kp is the cross section of a unitary scat- 
tered tTinei/fmax is the conditional probability of inelastic 
scattering if an electron hits the impurity, and 



1 

Thit 



2nj 

■KV 



(24) 



describes the typical 'hitting rate'. 

As in the WL experiment discussed above, fits to exper- 
imental data have to be done with the e-independent de- 
phasing rate. For a comparison with experiment we there- 
fore have to give the e-independent dephasing rate which 
for k = 1 is obtained by solving the equation 

LJT,B,L) k 



~ ST f de ~~^ < 



m L v (e,T,B) 



cosh 4 (e/2T) 



(25) 



where L v = \JDt v is the dephasing length. Notice that 
the actually measured dephasing rate depends on the length 
of the ring. It also differs from the WL result due to the 
different energy averages. 



IV. NUMERICAL RESULTS FOR DEPHASING RATES 
FROM AHARANOV-BOHM OSCILLATIONS 

In order to evaluate the dependence of the dephasing 
rate on magnetic field, temperature, and ring length from 
Eq.(27), we require the T-matrix for the single impurity 
Kondo model defined in Sec. II. By using the equation of 
motion method this can be expressed 



r v (e), (21) T a {u>,T,B) = J{S Z ) + J 2 (Sc^; Sowcv) (26) 

where ((...)) denotes a retarded correlation function and er 
are the Pauli spin matrices. We calculate Eq. (|26p by apply- 
ing the numerical renormalization group (NRG) method 31 
for finite temperature dynamics 32 . At finite magnetic field, 
it is also important to use the reduced density matri)^ 3 - to 
evaluate the above dynamical quantity. For all calculations 
presented here we used a discretization parameter for the 
conduction band of A = 1.5 and we retained 960 states per 
NRG iteration. We checked that this number of states was 
sufficient to maintain particle-hole symmetry of the spec- 
tral densities ImT t (w,T,B) = ImT;(-o;,T, B) at this rel- 
atively small value of the discretization parameter. The 
Friedel sum rule for the T = spectral density was satisfied 
to more than 1% accuracy in our calculations. 



Fig.|5]shows the numerical evaluation of y/ 'Ai(B) for var- 
ious T at a given length L = lOLhit where Lhit = V Diy^t ■ 
Notice that, if lengths are measured in units of Lhit, the 
amplitude (AiT^T^) 1 / 2 becomes a universal function of 
B/T K and T/T K . 

Fig. [TO] gives the magnetic field dependence of the de- 
phasing rate at various T and a fixed ring length L = lOihit- 
For large magnetic fields, B ^> T, Tk, the dephasing rate 
is expected^ to vanish proportional to (T/B) 2 / ln 4 [B/T K ], 
consistent with the numerical results (the precise form also 
depends on L/L v , see below). 

Fig.rjTJshows the results for 1/t^s as a function of T for 
different strengths of the magnetic field B at a fixed ring 
length L = lOLhit- While the maximal dephasing occurs 
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FIG. 9: Amplitude of the Aharonov-Bohm oscillations [in 
units of ((27r) 3 /(rhitT K )) 1/2 ] as a function of the applied mag- 
netic field, B (in units 7k, /iB/fe = 1), for different temper- 
atures T (in units 7k) obtained from NRG calculations for 
L/Lhi t = 10. The rapid rise of the oscillation amplitude re- 
sults from the suppression of dephasing (see Fig. 1 10 p by po- 
larizing the spins. Here, and in the remaining figures in this 
section, the symbols represent the discrete values of (B,T) at 
which NRG calculations were carried out. 
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FIG. 10: The dephasing rate as a function of the magnetic 
field, B (in units 7k) for various temperatures, T, and a ring 
of length L = lOLhit- Note the rapid suppression of 1/t v 
especially for low temperatures. 



for T — T K for small fields B < T K , it shifts to larger val- 
ues (T ~ B) for B ^> 7k. For high temperatures, small 
fields and not too large L/L v , see below, T ^> 7k, B the 
dephasing rate is well described by the Nagaoka-Suhl for- 
mula^ 1/t v {T) = 2^7 w 2 3/ 4 +1 3 n 4 4 T/TK - 

Fig. [T2l shows the dependence of l/T Vt s(T) on the ring 
length, L, at zero magnetic field. As pointed out above, 
the L-dependence of the experimentally measured dephas- 
ing rate enters through the energy average of Eq. (|25|) . As 
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FIG. 11: NRG results for the dephasing rate as a function of 
T and for different values of B. Thit = j^- as defined above, 
T and B are given in units of 7k^ (we set ^b/&b = 1). The 
values are for a ring of length, L, L = 10-Lhit where L^it = 
V Drmt- While the maximal dephasing occurs for T ~ 7k for 
small fields B <C 7k, it shifts to larger values (T ~ B) for 
B > 7k. 



can be seen from Fig. [12] the dephasing rate only changes 
by a factor 1/4 on increasing the ring length by a factor 
200. We included curves for the rather academic cases 
L = 100 — lOOOLhit (the amplitude is too strongly sup- 
pressed to be observed for such large ring lengths) in order 
to show this weak L-dependence. Fig. [12] also shows that 
for temperatures T < 7k the dephasing rate becomes en- 
tirely L-independent. This reflects the fact that the energy 
resolved dephasing rate, l/r^^e, T), establishes a deep 
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FIG. 12: Temperature dependent dephasing rate for various 
ring lengths L, measured in units L^it = \J ^r^- and calcu- 
lated via NRG for B — 0. The logarithmic suppression with 
L for larger T arises due to the interference of electrons with 
energies larger than T. 



minimum at e = for T < T< and the integral on the 
right hand side of Eq. (f25|) is therefore well approximated 
by setting e = 0. To be precise, in the limit of ring lengths 
L ^> L v the integral on the right hand side of Eq. (|25|) (for 
fixed T and B) is dominated by the saddle points of the 
function 



L 



-In 



4 4) ( £ ) 



4 In 



cosh 



V2T 



For temperatures T < max{S,TK} / has a saddle point 
at e = 0, which for temperatures T > Tk becomes unsta- 
ble. At very large ring diameters T/Thit 3> 10 2 , a second 
saddle point at e = ^ — ; — r T starts to dominate the 



2Lh it In 2 



integral for T > Tk. Here = \/Dt^ is the diffusive 
length scale corresponding to the time = and intro- 
duced above. Although this limit is rather academic it is 
interesting that for such big rings dephasing is dominated 
by rare events of highly excited thermal electrons scatter- 
ing from the magnetic impurities. This originates from the 
fact that high-energy electrons scatter less effectively from 
Kondo spins, as Kondo renormalization becomes less effec- 
tive for e ^> Tk. Inserting this second saddle point into 
Eq. ([25]) one finds that the length dependence of l/r„ for 
high temperatures follows 



(27) 



r v (T,L) 



hr 



2L hit In 2 



explaining the weak suppression of 1/t v for large ring 
lengths shown in Fig. [121 
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FIG. 13: Amplitude of the AB oscillations in units of e 2 /h as 
a function of T/T K for T = 40 mK (A) and 100 mK (o) mea- 
sured by Piere et al^£, assuming Tk = 10 mK. Solid and 
dot-dashed lines are the numerically calculated amplitudes 
with fitting parameters described in the main text. As for 
very high magnetic fields, B 3> 100 Tk, numerical errors in- 
crease when the dephasing rate becomes very small, we used 
an extrapolation of the numerical results, 1/t v oc l/B 2 , in 
this regime. The saturation of the amplitude at these high 
fields arises as the dephasing due to electron-electron inter- 
action dominates. The data is equally well described by the 
fits used in Refs. [2(| and [fj (dashed lines), see main text. 
For the solid lines we used the same values for the dephasing 
rates (r ee = 5.4 ns and 9.9 ns for T = 100 mK and T = 40 mK, 
respectively) as in Ref. [g], where r ee oc T -2 / 3 was assumed. 
For the dot-dashed curve we use instead r ee = 13.5 ns for 
T = 40 mK since one expects theoretically^ that r ee oc 1/T 
for L 3> Lff, (note, however, that L ~ L<j, in this regime ex- 
plaining the rather weak dependence on r ee ). 



V. DISCUSSION AND CONCLUSIONS 

In this paper we generalized previous results for the de- 
phasing rate due to diluted Kondo impurities as measured 
in the weak localization experiment to describe dephasing 
due to arbitrary diluted impurities. Furthermore, we investi- 
gated how magnetic fields modify the dephasing rate due to 
Kondo spins as can be measured in mesoscopic Aharonov- 
Bohm rings. We give results for the numerically evaluated 
dephasing rate as a function of the magnetic field, temper- 
ature, and the ring length. 

The influence of magnetic impurities on dephasing has 
been studied in a number of magneto-resistance experi- 
ments in Cu, Ag or Au wires doped with magnetic impuri- 
tiesi 2 ' 3 ' 5 ' 15 ' 16 More recently, high-precision experiments us- 
ing ion-implanted Fe impurities in Ag wires allowed a quan- 
titative comparison with our theory for spin-1/2 impurities, 
see Refs. [15] and [16] for a critical discussion. Such studies 
using samples doped with magnetic impurities have, to our 
knowledge, only been performed in the spin-glass regime 36 
using magnetic ions with tiny Kondo temperatures. Pierre 
et a/^S studied rings made from nominally clean Cu wires. 
As these wires show a saturation of the dephasing rate 



(determined from weak localization) at low temperatures 
(30mK< T < IK), it was suspected that tiny concen- 
trations of magnetic impurities with Kondo temperatures 
below 30 mK may be at the origin of the observed satu- 
ration. This picture could be confirmed as measurements 
of the amplitude of Aharonov-Bohm oscillation displayed a 
dramatic rise by almost an order of magnitude in moderate 
magnetic fields (see Fig. [T3|) . proving the magnetic origin 
of the \ow-B, low-T dephasing. 

As neither the concentrations nor the type(s) of magnetic 
impurities are known, a parameter-free comparison to our 
predictions is not possible for these systems. Assuming Mn 
impurities, believed to be characterized by a Kondo temper- 
ature of the order of 10 mK, 37 and, using the same dephas- 
ing rates due to electron-electron interactions as in Ref. [5] 
(r ee = 5.4 ns and 9.9 ns for T = 100 mK and T = 40 mK, 
respectively) we obtain the fits shown in Fig.[T3]for a g fac- 
tor of g » 1.4 and an impurity concentration of 2.7 ppm. 
We have also added a curve at T = 40 mK (dot-dashed line) 
which uses r ee = 13.5 ns (keeping all other parameters iden- 
tical) to take into account that one expects theoretically 35 
r ee oc 1/T. The fits and the extracted parameters are not 



very reliable as can be seen from the observation that the 
data has been equally well described in Refs. [5] and [20] by 
the simple perturbative formula ^^^a^ = ■ \ ^ B ^ k ?i T ^ 

r r T Vi g{B) smh(g fi B / k B T) 

with g = 1.08, see Fig. El 

For a more meaningful comparison to our results, exper- 
iments on AB rings, doped with magnetic impurities with 
a higher Kondo temperature would be highly desirable. On 
the one hand, such experiments and their theoretical inter- 
pretation can reveal basic dephasing mechanisms in metals, 
on the other hand, they can be used to obtain insight into 



the physics of strongly correlated dynamical impurities and 
their interactions. 
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